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A new method is introduced and used to calculate the statistics of the 
microscopic fluctuations of charge carriers in devices. By expressing the 
fluctuations of the carriers in terms of elementary transfer fluctuations, we 
are able to separate the induced fluctuations from the spontaneous fluctua- 
tions experienced by such carriers. This enables us to treat correlation 
effects in the dynamical portion of the problem and reserve for the statistical 
portion only well-defined, uncorrelated random forces whose statistics are 
readily calculated. The method includes all important correlation effects as 
well as multiple-decay-time relaxation effects and, thus, it fills a gap in the 
Langevin method as well as in the impedance-field method of calculating 
noise in devices. The method is suitable for treating non stationary as toell 
as stationary noise, and in some cases can be used directly on macroscopic 
problems. We also present a derivation of a recently introduced expression 
for diffusion noise of carriers -whose mobility is a nonlinear function of 
applied electric field. This microscopic approach may further illustrate 
the origin, nature, and treatment of fluctuations in devices. 

I. INTRODUCTION 

In this article we describe a very simple means of treating micro- 
scopic fluctuations in noise theory. The method is simple in the sense 
that it focuses attention directly on the heart of the matter, the 
elementary processes which give rise to device noise. It is also simple 
in the sense that no sophistication in probability theory is used beyond 
an understanding of simple shot noise. Nonetheless, the method is 
rigorous under the rather mild constraint that the fluctuations are 
sufficiently small that the equations governing the noise are linear. 
The method has the added advantage that it can be used nearly as 
easily for nonstationary noise as for stationary noise. We make no 
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claim that this method is an advance in the philosophy of noise ; we do 
claim, however, that it is adequate for solving many noise problems 
of practical interest. 

There are two equivalent 1,2 methods of calculating device noise, the 
Langevin 3 method (LM) and the impedance-field 4 method (IFM). 
Both methods are characterized by an inherent simplification : namely, 
the separation of the task of calculating the spontaneous fluctuations 
of the current carriers in each elemental region of the device, and the 
task of calculating the observable response to these fluctuations at the 
external contacts of the device. The former task, the treatment of the 
microscopic fluctuations, is simplified because in dealing with the 
source of the fluctuations one can focus attention on the statistics of 
the microscopic variations inherent to the local physical conditions, 
which in turn are determined by the (noiseless) state of the device 
during operation. As a result, both the LM and the IFM are primarily 
concerned with the latter task, the coupling of the microscopic fluc- 
tuations to the macroscopic, observable voltages and currents. This 
task is also well-defined because the influence of the carriers in one 
region of the device on the carriers in another region, and on the 
contacts, has often been studied in detail in attempts to understand 
the dc and ac operation of the device. Thus, it is important to comple- 
ment the LM or the IFM with the microscopic method described 
below. When this is done, it can be claimed that in most cases, if one 
understands the device sufficiently to calculate its noiseless operation, 
one can calculate the device noise as well. This should be of assistance 
to those tackling the noise in new and/or unfamiliar devices from 
scratch. 

For example, if a device is sufficiently well-understood to be charac- 
terizable by an equivalent circuit, one can often introduce equivalent, 
random voltage and current sources to simulate the noise in the 
device. 6-8 Using such sources, a circuit designer with little interest in 
noise theory can readily calculate the size of the noise in the circuit 
employing the devices of interest. In a similar way, a person working 
with individual devices may find it convenient to have a simple scheme 
to translate the physical processes with which he is familiar into noise 
sources and to quickly evaluate their effect on the performance of the 
device of interest. 9-10 It is hoped that the method presented here will 
be used in such situations. 

It is important to realize that charge carriers in any device fluctuate 
in response to random forces exerted on them. 3 The response to a 
specific impulse continues in general long after the impulse causing it 
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has ceased, and, in the interim, subsequent impulses will further alter 
the induced fluctuation. In addition, a fluctuation in the distribution 
of one type of carrier can induce fluctuations in the distributions of 
other types of carriers. Thus, whereas the statistics of the spontaneous, 
random forces may be quite simple, those of the carriers can be some- 
what complicated owing to the correlation between the various in- 
duced fluctuations. The key to the simplicity of the method presented 
here results from the separation of the correlation effects from the 
statistical problem. These correlation effects are not neglected. Rather, 
they are included in the dynamical portion instead of in the statistical 
portion of the treatment. As it turns out, this results in the primary 
simplification achieved with our method. 

In what follows, we shall use several examples to introduce and 
elaborate our microscopic approach. The first example, the decay of 
charge stored on a leaky capacitor, will motivate the method and 
illustrate how this technique can be used to treat certain macroscopic 
problems as well as microscopic ones. The second example concerns 
transfers between a two-level system. Here correlations are of primary 
importance, and our method is seen to treat these adequately yet 
simply. The third example illustrates how velocity fluctuations can be 
decomposed into transfer fluctuations, which are much easier to treat. 
Complicated scattering mechanisms including multiple decay times, 
which are not normally covered in the usual Langevin 3 or impedance- 
field 4 methods, can be handled with relative ease. The fourth example 
considers recombination to illustrate how correlation effects can be 
treated efficiently and effectively. By working in the time domain, we 
can see how the method works for nonstationary 9 as well as for station- 
ary statistics. The statistics are treated in detail in the appendices. 
In particular, a recently used expression 11 for the diffusion noise of 
carriers having a nonlinear mobility is derived. 

It is hoped that our discussion will provide, for the nonspecialist, 
further insight into the physical nature and mathematical representa- 
tion of noise in general. 

II. SOME SIMPLE EXAMPLES 

Since the primary purpose of this paper is to elucidate a general, 
practical approach to solving the microscopic portion of noise problems, 
it seems best at first to describe this approach in terms of simple 
examples. Although this method is best suited to treat noise at the 
microscopic level, for purposes of illustration we shall commence with 
a macroscopic example. 
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Fig. 1 — (a) RC circuit in which charge initially placed on the capacitor decays to 
ground through the resistor, (b) Noise equivalent circuit of (a). 



Consider the RC circuit shown in Fig. la in which we assume that 
capacitance C is a function C(Q) of the charge Q, which it stores. Let 
us assume that initially (t = 0) a charge of size Q° is stored on the 
capacitor. Subsequently (t > 0) the charge will decay away through 
the resistor R. In the absence of noise, the charge as a function of 
time Q°(t) satisfies the equation 



dQ°(t) 
dt 



Q°(0 



RCLQ°(t)T 



(1) 



subject to the boundary condition that Q°(0) = Q . 

We know, however, that the thermal (Brownian) motion of the 
electrons in the resistor gives rise to a noise current. Thus, rather than 
the uniform charge decay predicted by (1) for noiseless conditions, the 
charge decay is in fact somewhat random. Moreover, the noise is not 
solely governed by the resistor. If, during a time interval, the noise 
current in the resistor is such as to draw too many charges from the 
capacitor, then the voltage on the capacitor will decrease and the 
subsequent current will be reduced. This is just the response of the 
circuit to the noise generated by the resistor. In such a problem, we 
would ordinarily determine the statistics of the noise current i(t) 
associated with the resistance itself, and then determine the fluctua- 
tions in the charge Q(t) on the capacitor in terms of i(t). We would be 
implicitly assuming, and rightly so, that only the resistor, and not the 
circuit exterior to the resistor, determines the statistics of i{t), the 
thermal noise current generated by the resistor. We would then relate 
Q(i) to i(t) using an equation of the form 



dQ 

dt 



= -i(Q) - m, i(Q) = 



Q 



RC(Q) 



(2) 



(The noise equivalent circuit is shown in Fig. lb.) We would then 
write Q{t) as a sum of its stochastic average Q°(t) and a fluctuation 
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8Q(t) = Q(t) - Q°{t); assume that (8Q{t) 2 )l « AQ, where 

AQ = (dI/dQ)/(d*I/dQ*); 

and expand eq. (2) to first order in 8Q(t), obtaining a nonlinear equation 
[eq. (1)] for the Q"(t) and a linear equation, 

d8Q 8Q 

~W=~W)~ t(0, (3) 



for 8Q in terms of i(t). In eq. (3) r(t) is defined by 

(4) 



r(0 MflC«2)/l<2=QV> 

Solving (3) for 8Q(t) in terms of i(0 permits the statistics of 8Q(t) to 
be obtained from those of i(t), which for now we assume we know. In 
this way, we can determine the noise associated with the variable of 
interest Q(t) from a knowledge of the noise associated with a more 
simply characterized noise variable i(t). 

The foregoing procedure has one drawback : one must be extremely 
careful, in general, to properly introduce such noise terms as i(t) into 
otherwise noiseless dynamical relations. As we shall see in the examples 
considered below, there is a one-to-one correspondence between the 
transfer processes characterizing the problem of interest and the noise 
terms which one introduces. Thus, if one writes down several coupled- 
rate equations, each involving several transfer processes, and introduces 
but one noise term per equation, one finds, in general, that these noise 
terms are not simple, being correlated statistically to one another. 
Since, for simplicity, one would desire noise terms to be uncorrelated, 
some care must be used in including them in the rate equations. We 
now outline a procedure in the context of the above example that can 
be used in more complicated problems to insure that such noise terms 
are included properly. 

The decay of charge from capacitor to ground is accomplished by 
transfer of individual electrons. Let <, be the time at which the ith. 
electron leaves the capacitor. Then we may write the following kinema- 
tic relation for the decay of the charge on the capacitor, 

^f=- 9 D«-(,), (5) 

where q is the size of an elementary charge. If we consider an ensemble 
of RC circuits of the type shown in Fig. la, in which an initial charge 
of Ql is decaying, each decay will be characterized by a different set of 
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times [ti\ at which elementary charges leave the capacitors. Each U, 
therefore, is a random variable whose probability distribution is, in 
general, dependent upon all l, for the preceding events (tj < U). A 
completely rigorous derivation of the noise would, therefore, involve 
calculating the probability of each such sequence of times [tt], 
including all correlation effects, and then using these probabilities to 
ascertain the statistics of the noise. Were it not for these correlation 
effects, the ti would be independent, and the statistical problem would 
be greatly simplified. 

We shall now recast eq. (5) into a form that greatly simplifies the 
correlation problem by decomposing the current into a spontaneous 
portion d(t) and an induced portion R(t). The spontaneous portion is 
governed by the sources of the noise, and the induced portion is con- 
trolled by the instantaneous state of the device, in this case the stored 
charge Q(t). Returning to the charge decay problem, we rewrite eq. (5) 
in the following form : 

^ = -qR(t) - qd(t), (6) 

where 

d(t) = Z<Ht - U) - R(t), (7) 

and where R(t) is the "dynamical" rate of charge loss. By "dynamical" 
rate we mean that R(t) is, in general, a function of (that is, is deter- 
mined by) the dynamical variables of the problem. In this case, we 
have 

If there is a fluctuation 8Q in Q, then a fluctuation 8R(t) occurs in 
R (t) also, which in this case is given to first order in 8Q (t) : 

q8R(t) = 8Q(t)/r(t), (9) 

where t is defined in eq. (4). Returning to eq. (6), if we write 
Q — Q° + 8Q as we did in passing from eq. (2) to eq. (3), then we 
obtain an equation for the noise 8Q [which corresponds to eq. (3) 
above], 

d8Q 8Q 

-Jf = " W) " qd(t) - (10) 

In this case, qd(t) corresponds to the noise current i{t) generated by the 
resistor, which, as discussed above, is independent of the 8Q associated 
with the capacitor. Thus, d(t) serves as a statistical driving term. We 
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may calculate 8Q(t) in terms of d(t) and from the statistics of d(t) 
determine those of 8Q(t). 

Let us return to eq. (7) for a moment to point out why d(t) is to be 
regarded as the spontaneous portion of the current. The — e £«' 8 (t — t % ) 
term is, of course, the entire current associated with the charge decay. 
The dynamical rate term Rit), however, is a function only of the dy- 
namical variables of the problem and does not contain noise sources. 
In our example, R(t) involves only the charge Q(t) = Q°(t) + 8Q(t), 
where 8Q(t) is the charge fluctuation induced by the noise sources 
acting on the device. R(t) does not involve the noise sources them- 
selves. Thus, if d(t) = 0, 8Q, the response, vanishes. Hence, inasmuch 
as d(t) is the difference between the total current and the noiseless- 
plus-induced portion of the current, it follows that d(t) can contain 
only the spontaneous portion of the current. The advantage of starting 
with eq. (5) and proceeding as we did to eq. (10), rather than attempt- 
ing to write eq. (2) [or eq. (3)] a priori, will become evident when more 
than one process is involved in the problem. By introducing a dt(t) 
and Bi(t) for each process i, and noting that each dt(t) involves only 
spontaneous fluctuations and, hence, must be independent of all the 
other dj f j 5^ i, we can readily express the fluctuations of interest in 
terms of the independent statistical driving terms d». 

If we solve eq. (10) for 8Q(t), we find that 



8Q 



(t) = £ dl'exp [- J'dt"/r(t")] L-qd(t'n (11) 

Often one is most interested in the mean-square fluctuation (8Q 2 ) for 
some time t. From (11) this is given by 

W(0> - 9 2 /' dt[ f* dfc'exp [- f* t dt"/ T (t")] 

X exp I" - j' t dt"/r(t")] {d{t[)d(t' 2 )). (12) 

If we know the statistics of d(t) (in this case, those of i(t)/q associated 
with the resistor), we can calculate (d(ti)d(tt)), and hence (8Q 2 (t)). 
Note that, in general, we must also know the noiseless solution Q°(t) 
[eq. (1), and also see eq. (4)]. For simplicity here, let us assume that 
i(t) is pure thermal noise so that 

5 2 (d(«i)d(< 2 )) = 2kTG8(h - U), G = l/R. 

It follows then that 

(8Q»(0> = kTC, (13) 
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the usual result. [Note: (6Q (h)8Q (t 2 )) = kTC exp (- \h - t 2 \/r). This 
illustrates how eq. (10) maintains the correlation between fluctuations 
at two different times, while the fluctuations themselves are driven by 
a source without correlation, 8(t — t')r\ When we turn to purely micro- 
scopic processes, we shall find the statistics of d(t) are governed by the 
rate function R(t). 

In what follows, the very important distinction between spontaneous 
and induced fluctuations will be used repeatedly. The advantage of 
this macroscopic example is that the separation between the two may 
be clearly visualized. From the point of view of the RC circuit, fluctua- 
tions in the current generator are spontaneous and induce fluctuations 
in the charge decay, 8Q. In addition, the rate term R(t) is a function 
only of Q and clearly includes only the induced and not the spontaneous 
fluctuations. The reader may find it helpful in subsequent examples to 
refer back to this simple model to clarify the somewhat more subtle 
distinctions between induced and spontaneous fluctuations at the 
microscopic level. 

We now consider another simple example, this time a truly micro- 
scopic one. Let us consider two states a and b containing n„(t) and 
nb{t) charges, respectively. [These states, for example, might be two 
regions of phase space (dxdv) :' one region for (x a , v a ) and one for 
(x&, v 6 ), or a might be trapped electrons and b free electrons.] Let us 
assume that charges are flowing to b from a at a rate Rba = Rba(n a , n&) 
and to a from b at a rate R a b = R a b(nb, n a ). If charges leave a and enter 
6 at times hai and leave b and enter a at times t a u, then by analogy with 
(5) we write 

n a = -Z*(t- had + £«(«- t abi ) (14a) 

t j 

and 

n b = - L 8(t - t abi ) + ZHt~ had- (14b) 

t t 

Following the previous example, we rewrite (14) in the form 

n a = -R ba (t) + Rab(t) - d ba (t) + d ab (t) (15a) 

and 

n b = -Rab(t) + R b a(t) - d ab {t) + d ba (t), (15b) 

where, of course, 

d ab (t) = Z8(t- tabi) - Rab(t) (16a) 



* Note that boldface capital letters denote matrices ; boldface lower-case letters 
denote vectors. 
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and 

d ba (t) = Z B(t - t bai ) - Rba(t). (16b) 

If we write n a = n° a + 8n a , n b = n° b + 8n b , insert into (15), and expand, 
we obtain the following equations for the noiseless quantities n%, n° b : 

n°a = -Wait) + R°at(t) (17a) 

and 

nl = -R° ab (t) + R° ba (t), (17b) 

where 

R° ba = R ba {< {t) , n% (t) ), R° ab = R ab (n° b (t) , n° a (t) ) . 

For the noise 8n a , 8n b) we obtain the linear relations 

. . / 8R ba 8R ab \ s ( 8R ba 8R ab \ s 

Sn tt = — I -z ^ — 1 dn a - I -r r — I 8n b 

\ 8n a 8n a / \ 8n b 8n b / 

- d ba (t) + d ab (t) (18a) 



and 

( 8R a b 
5n, 



/ on a b 8R ba \ 1 8R ab 8R ba \ 

8n b = — ( -^ r — J 8n a — [ — r — I 8n b 

8n a / \ 8n b 8n a ) 



- d ab {t) + d ba (t), (18b) 



from which 8n„ and 8n b can be determined in terms of d ab and d ba . 
(The linear operators of the form 8R/8n are evaluated at their noiseless 
values.) Since there are only two states, and since n a + n b = constant 
(n a + n b = 0), it comes as no surprise that 8n a (t) = —8n b (t); one 
state's loss is the other's gain. Nonetheless, the source or driving terms, 
d ab and d ba , are independent, and the correlations between 8n a and 8n b 
are included in (18) through the presence of the dynamical terms. Since 
this is an important point, we shall discuss it more fully below. 

If we ignore for the moment the random nature of the flow of charges 
from a to b (and from b to a), then eq. (17) tells us that we have a 
"smooth" continuous flow of charges from a to b at a rate of R° ba {t) 
and from 6 to a at a rate of R° ab (t) . Noise enters the problem when we 
note (as we have above) that charges are actually transferred at times 
\t ba i) and {t ab i), the (ensemble) average rate of occurrence of these 
times being R ba (t) and R° ab (t) . Since R° ba (t) and R° ab (t) depend only on the 
steady-state solution (n° a , n° b ) to eq. (17), these quantities are not 
affected by the details of a particular set of fluctuations. This means 
that if R° ba and R° ab govern the statistics of the transfers from a to b and 
b to a, then the individual transfer times fd-}, [t ab i\ associated with 
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Rla and R° ab must be statistically independent. [If they were cor- 
related, then R° ba and R° ab would reflect this correlation much as do Rba 
and R a b in eq. (15). J In discussing the discharge of the capacitor, we 
pointed out how d{t) acted as a noise source containing only the 
spontaneous fluctuations. Above, we have noted that these fluctuations 
are governed by the average rates R°. This suggests that to calculate 
the statistics of d ab and d ba , which according to (18) are needed to 
calculate the statistics of 8n a and 8n b> we may write the d(t)'s in the 
following form : 

d ab (t) = £ 8(t - t° abi ) - R° ab (t) (19a) 

i 

and 

dba(t) = £«(<- Gat) - R° b a(t). (19b) 

i 

We stress at this point that for a specific event ({t ab i, t ba i\) we need not 
demand that the right-hand sides of (16) and (19) be equal. This is 
because we eventually average the dependence of n a and n b on d ab and 
d ba over all events in calculating correlation functions and spectral 
densities. All that is necessary is that the statistical properties of the 
two forms of d ab and d ba be the same, at least to lowest order in the size 
of the fluctuations. We shall see below that this is indeed the case. We 
are able to calculate simply the statistics of d ab and d ba because in 
writing (19) we have cast these random variables into a form in which 
the distinction between spontaneous and induced fluctuation no longer 
enters. We can do this because d involves only spontaneous (hence 
independent) fluctuations and, therefore, can be written in terms of 
independent events. 

With independent transfer times, it is straightforward to calculate 
the statistical distributions of d ab and d ba as given in (19). This calcula- 
tion is carried out in Appendix A. The autocorrelation functions for 
d ab and d ba , which are the spontaneous-fluctuation, noise-source terms, 
are given by 

(d ab (ti)d ab (t 2 )) = R° ab (ti)Hh ~ U) (20a) 

and 

(d ba (h)d ba (t 2 )) = RKtiMh - h), (20b) 

as expected for pure shot noise associated with independent events. 
In addition, if R° ab and R° bn are independent of time (so that the noise is 
stationary), then the spectral densities of d ab and d ba are given by 

SUf) = 2/C (20c) 
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and 

SUf) = 2iC. (20d) 

The form given d a b and db a in eq. (19), or rather their statistical form, 
eq. (20), is the basic assumption that enters our approach. That it is 
valid so long as the fluctuations 8R in the rates R = R° + 8R satisfy 

(5# 2 >* « R° (21) 

may be motivated by the following. If a fluctuation Bn in n leads to a 
fluctuation 8R in R, then the statistics of the immediately following U 
will be governed by the altered rate. This suggests as a first-order 
iteration in the fluctuation that 

(d(h)d(ti)) « R(h)S(ti - h). (22) 

However, R = R° + 8R ; thus, the error in (20) is of the order of 
8R/R", which we assume to be small. A similar argument applies to 
residual correlation effects among the d's. Thus, so long as (21) is 
satisfied, we expect (20) to be valid to order 8R/R". Since we are 
discussing a linear theory, this is good enough for our purposes. 

Returning to (18), assuming stationary noise so that we may use 
Fourier-transform techniques, and representing that the linear opera- 
tors are simple decay rates for purposes of illustration, we may write 

iu8n a = —8?i a /T a + 8)ib/Tb — db a + d„b (23a) 

and 

iu8?ib = -\-8n a /T a — 8nb/Tb — (Lb + dba, (23b) 

which become, upon solving for 8n a (u) and 8nb(o}) in terms of db a (o>) 
and daft(w), 

Bn„ = -dm = (d ab - d ha )(iu + l/r)' 1 , (23c) 

where 

l/r = l/r. + 1/n. (23d) 

Since the d a b and db a are independent, for spectral densities we obtain 

S s a n (f) = Sl n (f) = | fa + l/r\- 2 lSUf) + SUf)2 (24a) 

= \iw + 1/t|- 2 (2JC + 2R° b a) (24b) 

= SM) = -Sfim (24c) 

where the expressions in (24c) are cross-spoctral densities (/ = co/27r). 
In such a problem as this with transfers permitted only between two 
states, the correlation between 8n a and 8nb is maximal. (Introducing 
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more states will lessen the correlation between any two.) Nonetheless, 
the important features of the correlation are clearly obtained in the 
result (24b) : the small-signal decay rate (1/r) is the sum of the separate 
decay rates, and the spectral density is proportional to a full shot- 
noise term derived from the sum of the noiseless rates (Rab + Rio). 
This example is admittedly very trivial and could have been done 
more easily by other means. In the following sections, the advantages 
of the present method will become evident. Note, however, that there 
was no question as to how to include the driving terms when writing 
down the coupled noise equations (18) and no question regarding the 
independence of these driving terms when it came to calculating 
spectral densities in (24). 

Let us pause before going on to summarize the logic we have used in 
arriving at our method. First we recognized that noise arises owing to 
the randomness in the times at which charge carriers change their 
state. Then for each transfer process, we separated the expression for 
the particle current flowing between any two states into a dynamical- 
rate term R(t) that includes all induced fluctuations, and a driving 
term d(t) that includes only spontaneous (and therefore independent) 
fluctuations. In this way, we are able to shift the statistical correlations, 
so to speak, out of the statistical term and into the dynamical term. To 
calculate the statistics of our driving term, we noted that since the 
fluctuations were characterized by a mean rate of R°(t), which depended 
only on the noiseless solution, we could rewrite d(t) in a natural way 
in terms of independent events. [These latter events cannot be coupled 
since their statistics are governed only by R°(t) and not by preceding 
events.] Although this step represented an approximation, we argued 
that it should be all right for small fluctuations. With the new form 
for d(t), its statistics were readily calculated. Finally, knowing the 
statistics of d and the linear relation between Sn and d, the statistics of 
8n including all correlation effects of interest, which we desire, are 
straightforward to obtain. In the next two sections, we shall carry out 
the above procedure for two examples from start to finish. 

In passing, we make reference to Lax's discussion 12 in which he 
showed the basic limitation of any source theory of noise. The heart of 
Lax's argument is that one cannot expect to model on the time scale 
of the duration of a spontaneous fluctuation (during which time the 
system appears to be reversible) by using the essentially irreversible 
dynamics embodied in eq. (18). As for time scales on the order of the 
response time of the carriers to spontaneous fluctuations, however, all 
is well. Thus, as Langevin probably recognized long ago, the price paid 
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for using source terms is small. The whole trick is to include them 
unambiguously. 

III. NOISE IN CHARGE TRANSFER 

The purpose of this section is to stress the insight one can gain by 
starting with the most elementary transfer processes taking place in a 
device. In this case, we avoid making assumptions about the statistics 
of the velocity fluctuations of carriers by calculating them from oc- 
cupancy fluctuations, which can be understood much more simply 
using our method. For simplicity, we shall ignore correlation effects, 
which in fact are shown to be negligible for the problem we discuss. In 
the next section, which is on recombination, we shall stress correlation 
effects and how they can be dealt with using the microscopic method. 

One of the most basic approaches to the problem of the storage and 
transport of carriers in a device is first to assign a density function to 
each carrier state. For example, n a (x, v, t)dxdv represents the number 
of carriers of type a at x with velocity v in the element dxdv of phase 
space. Similarly, n b (x, t)dx represents the number of carriers of type b 
trapped at x in volume dx. In this example, let us ignore trapping and 
recombination effects and focus attention on the transport via scatter- 
ing of a single type of carrier. This will lead to a general expression for 
diffusion noise which we can compare with the results obtained by 
Langevin and by Shockley et al. using the Langevin method and the 
impedance-field methods, respectively. 

We proceed as follows. If at times Uti a particle is scattered from 
state j to state i, then we have (by analogy with our treatments above) 
that 

**. -ZSit-tj.d + ZHt-t,,:) (25) 

at ji ji 

and 

tt - - E fi*(0 + E ««(*) - £ W) + £ <M0- ( 26 > 

at j j j j 

The subscripts *, j each designate a particular region of phase space 

dxdv at (x, v). (We do this to simplify the notation.) In (26), d ti (t) is 

defined by 

da(t) =■£*(<- Un) - Rait). (27) 

i 

It is important to note that we introduce a source term for each 
elementary process, that is, for each transfer process indexed by the 
ordered pair (i, j) . 
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As before, we let n,- = w° + 8tii, insert into (26), and expand to 
lowest order in 8nt. The noiseless equation that results in 

^T= -X.R°n(t) + ^R%{t) (28) 

is just the Boltzmann equation, R°j being a function of n° and n°. 
Using x, v notation, (28) becomes 

dn(x, v, t) _ dn(x, v, t) dn(x, V, <) ■ F d?i(x, v, I) 

d£ dt dx m dv 

= - /dx' dv'(R(x, v; x', v') - 22 (x', v'; x, v)), 

the more usual form of the Boltzmann equation, where 22(1, 2) is the 
average rate at which particles are scattered from 2 to 1. In general, 
12(1, 2) is a function of ni and n 2 . Usually 22(1, 2) is taken to be propor- 
tional to ri2, however. To determine such noiseless quantities as 
mobility, etc., (28) must be solved for the noiseless n°(x, v, t). 
The equation for the noise 8n(x, v, t) = 8ni obtained from (26) is 

- Z «*) + L du(t). (29) 

i y 

In general, this is a coupled (linear, integral-differential) equation 
among the 8n,-. [JAs before, the linear operator 8R(n)/8n is evaluated 
at its noiseless value by inserting the noiseless solutions n° of (28).] 
For our purposes here, we simplify (29) by assuming (i) that jR>, is a 
function of (and not an operator on) the iik, and (ii) that in (29) we can 
ignore 8n k where k j& i. Assumption (i) is in fact the usual situation 
one has with the Boltzmann equation. Assumption (ii) can be made 
plausible in the following manner. By ignoring in (29) 8n k , k 9± i, we 
are ignoring correlations among the fluctuations 8rii and 8n k . Let us 
suppose that the scattering rates between i and k are equal for all 
N{N -j- 1) (i, k) pairs. Then a fluctuation 5n, will lead to a 8n k on the 
order of 8rii/N. The effect of 8n k back on state i will be of the order 
8n k /N . Thus, summing over the ./V states k ?± i, the correlated contribu- 
tion to 8rii which we lose by ignoring the 8n k in (29) is of the order of 
8rii/N [_ = N-(8rii/N)/N~], which for large N is entirely negligible. 
While we have assumed that the transfer rates between all i and k 
states are equal, a similar argument applies so long as the scattering 
to a few states is not favored. In this case, the 8n k dependence of these 
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few states must be included in (29). (We have already considered the 
two-state problem where correlation is largest. In our next example, we 
shall consider recombination in which correlation among several states 
is important.) Finally, we should note that our approximation is not 
equivalent to a relaxation-time approximation, which might be used to 
simplify (28). 

Returning to our example, and based upon the reasoning given above, 
we approximate (29) by 

Sfii = — E -Tzr 8n > + E t~^ 8n < 

j on, j dm 



and 
where 



and 



- E dM + E da(t) (30a) 

.i j 

8n { = -6n s /Ti(t) - d t {t) + d\{t), (30b) 

^^(t-tO* (3oc) 

d°i{t) m E *(/ - t»i) - E RiS), (31a) 



dm = EHt- tm) - E Rvtf). (3ib) 



The superscripts on d { , "o" and "e," designate "out" and "in" scatter- 
ing, respectively. We are able to lump the driving terms together 
in this way because (30a, b) involve only 5/?, and because the residual 
correlation between d° and d'i is of order 1/N. 

To proceed further, we approximate (31a, b) in the manner discussed 
above [see (19)] and determine the statistics of dn,-. Using the general 
results derived in Appendix A for d and then in Appendix B for 2ra as a 
function of d, we find for the correlation function of 5n„ Snj the 
expression 

(Sn.UJdnjUo)) - 8i£R° t (m) + «Kout)](r,/2)e-i'i-'2i'% (32) 

where we have assumed that the noiseless rates R"(in) = E/^y an d 
R°(out) = J^jR'ji [and therefore r, defined in (30)] are independent of 
time. One further simplification is often justified. Under the stationary 
conditions taken here to arrive at (32), R°(m) = R° (out). Furthermore, 
if Rji depends only on //, and is directly proportional thereto, it follows 
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that 

R°t(out)Ti = n° t . (33) 

Thus, we find the expected result that 

(«ni(<i)Sra y (* 2 )) = 8g»&rl«r-ni/r,. (34) 

We shall now use this result to determine the statistics of the velocity 
fluctuations of the individual carriers. It is velocity fluctuations rather 
than occupancy fluctuations that are usually taken as basic. We shall 
see, however, that occupancy fluctuations are the more fundamental 
of the two. 

Diffusion noise arises from the velocity fluctuations that a charge 
carrier undergoes owing to its interaction with the material in which it 
is confined as well as with the other carriers. Since the current in a 
device is made up of the linear superposition of the currents carried 
by each carrier, and since each contribution is proportional to the 
velocity of the individual carrier, a calculation of the current correla- 
tion function (to determine device noise) requires knowledge of the 
velocity autocorrelation function for each carrier. This latter quantity 
we can calculate at once from (34). We proceed in the following manner. 

By definition, the velocity of a carrier is given by 

ld 3 v'v'n{x,v',t) 
v(x,fl =^- (35) 

jrfVft(x ) v',0 

If we recall that n = n" -f- Sn, we can write (35) in the form v = v° + 5v 
to determine that 

fd 3 v'v'n°(x, v', t) 

V(«»,i) = — (36) 

d 3 v'n°(x, v', t) 

and that 

fd 3 v'v'8n(x, v', t) 
8v(x, = *—- - , (37) 

\d 3 v'n°(x,v',t) 

where we have assumed that the total number of carriers within a 
volume element dx remains at the noiseless value. We make this as- 
sumption only because we are interested here in diffusion noise only 
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and not noise from fluctuations in particle density. It follows at once 
from (34) (assuming stationarity) and (37) that the velocity correla- 
tion function of interest is given by 

(8v,(x h t)8vj(x 2 , h)) 

/ il 3 v'v'iVjn°(xi, v') 

- 5(X, - X 2 )— J — —— - e -|' I -'2l/r(x,.V). (3 8 ) 

w°(xi) MVn°(xi, v') 
In (38), i and j now correspond to the components of the velocity 

V = (V x , Vy, V t ). 

When one is dealing with stationary noise, it is more convenient to 
work with spectral densities than autocorrelation functions. Using the 
standard definitions of the spectral density function, 13 one finds using 
(34) that 

S Sn (Xl, Vi, x 2 , v 2 , /) 

= 5(xi - x 2 )5(v x - v 2 )4 / e- llT( - Il - vl) n°(x h v x ) cos utdt 

Jo 

- afv *- \*(v v ^ ^(x^viMx^vO . . 

- «(xi - x 2 )5 ( Vl - v 2 ) -___^_j- (39) 

It follows that the spectral density of the velocity fluctuations is given 
by 

q (t -r f\ - &(* v^ /• f/ „ 4»°(x 1 , v)t(X!,v) VjVj 

S^(x u x 2 , /) - 5( Xl - x 2 ) J dv 1+a)V(xi>v) ^^y 2 , (40) 

an expression which could have been obtained directly from (38). 

The spectral density given in (40) has several interesting features 
that we shall touch on briefly. Ordinarily for the frequencies of interest 
in devices, cj 2 t 2 « 1, owing to the very rapid carrier-scattering rates. 
In this case (40) reduces to 414 

S Sv .ij = 8(x, - x 2 )4D, ; //t°(x 1 ), (41) 

where Dy, is the ij component of the diffusion tensor defined by 

D«f-g£<l*(«!)ai*(«i + 0>. (42) 

This expression (42) can be derived from (38) by noting that 

8x(t) = j dt'8v(t'). 

J —CO 
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The diagonal components of D, D {j = Z), satisfy 
/ dvv\n°(x, v)t(x, v) 



Di = 



I dvv\n°(x, v) 



fdvvh°(x, v)/n°(x). (43) 



In this form, we recognize that the first factor is mm(x)/q, while for a 
thermal distribution, the second factor is kT/m. Here, m(x) is the 
carrier mobility in the ith direction at x derived using a standard 
Boltzmann equation approach. Inserting into (41), we find the usual 
result for thermal noise, namely that 

S Sw (x h x 2 , /) = «(xx - r 8 )4fcZV(x)/fon'(xi)). (44a) 

In Appendix C we generalize (44) to the case of field-dependent 
mobilities n{E). The result when an effective temperature T can be 
defined is that 

*.<*, x„ /) = ife - x 2) 4 ( *£)SeOgB/*K>. ( 44b) 

We are now in a position to compare the results obtained here with 
those obtained by Langevin 3 and by Shockley et al. 4 using additional 
assumptions. For simplicity, we shall work in one dimension and ignore 
spatial variations. In the Langevin 3 method, one begins by decomposing 
the force acting on each carrier into two parts, a damping force propor- 
tional to the velocity and a stochastic force of zero mean. The purpose 
of the latter is, of course, to produce the random fluctuations in the 
velocity. Thus, one writes 

8i> = -8v/t + h(t), (45) 

where r = mn/q. One then calculates the spectral density of 8v in 
terms of that of h, obtaining 

S iv (f) = j^fi£. (46) 

If one assumes that Sh(u) = Sh(0) [white noise corresponding to 
totally uncorrelated h(t)"2, then since 

^ = <*•) = / dfS iv (f) = & (0)1/4, (47) 
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one finds that S h (u) - 4kT/mr. It follows then from (46) that 

_ 4kTr/m _ 4kTn/e . ) 



/ dtm°(v)y 2 - 



1 + C0 2 T 



2„2 



(48b) 



This result agrees with (44) for ut « 1. If, however, we are interested 
in a?, for which tar W 1, then this result differs markedly from (40), 
which we write as 

unless r(y) = r, a constant. This, however, is seldom the case in 
realistic scattering problems. We may rescue the Langevin approach 
if, in place of (45), we pass to the frequency domain to write 

iu)8v = —8v/r u + h u , (50) 

in which the effective damping 1/t m is frequency dependent; and now 

&.</> = TT^k- (51) 

Again, we choose S h (f) = S h (0) and note that as / -► we must obtain 
(48). Thus, S h (f) = 4kT/niT, where t = r(w - 0). Inserting into (51) 
we obtain 

.._ 4 ldvn°(v)v 2 /T 

VcT/mr J (52) 

If (52) is equated to (49), one can solve for l/r«, that is, for the ap- 
propriate damping term to be used in the Langevin equation. This 
illustrates a major defect of the Langevin 3 approach, quite apart from 
inserting h u in a rather arbitrary manner. If the appropriate t u is not 
known a priori, one must return to a more fundamental approach such 
as that given above. This is important, for if x u = n u E u , then 

Mu = (e/m){ico + l/rj)- 1 = (e/kT) I dvv 2 n(v)t™ + 1/t(»)]~ 1 . 

Setting t u = ri does not make (52) equal (49). Thus, the original 
Langevin idea is not internally consistent in general. 
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In developing the impedance-field 4 method, Shockley et al. focus 
attention on the trajectory of a given carrier and on the carrier's 
deviations from its noiseless trajectory. Nonetheless, when it comes to 
writing an expression for the autocorrelation function of the velocity, 
they must postulate an expression equivalent to (38). [See eq. (48) of 
Ref. (4).]] In more complicated problems, it seems best to have a 
fundamental approach from which such microscopic correlation func- 
tions can be derived. For example, were it necessary to include 
correlation between certain pairs of states, several characteristic decay 
times would be present in the correlation expression. These can be 
included in a natural way if we derive 5v from eq. (37), as we have done 
here. Using less fundamental approaches, one must often rely on 
analogy and intuition. 

IV. RECOMBINATION-GENERATION NOISE 

Up to this point, we have been concerned with deriving the statistical 
distributions of microscopic variables. From a knowledge of the de- 
pendence of the macroscopic quantities of interest on these microscopic 
variables, the statistical distributions of the former can be calculated 
from those of the latter. In some cases, e.g., thermal noise in resistors, 
the macroscopic current or voltage spectral densities can be obtained 
very simply using thermodynamic arguments without recourse to 
microscopic methods. We now turn to an example, that of recombina- 
tion noise, for which the power of treating correlations among fluctua- 
tions using the microscopic method developed in this paper becomes 
apparent. 

Let us consider the following model of recombination. Let n be the 
number of mobile electrons, p the number of holes, n t0 the number of 
unoccupied traps, n, p the number of traps containing a trapped hole, 
and n ln the number of traps containing trapped electron, all within a 
unit volume. Recombination occurs when a trap containing an electron 
captures a hole or when a trap containing a hole captures an electron. 
In Fig. 2, we present schematically the eight different processes which 
enter into this recombination-generation model. In Fig. 2, as in the 
equations for the dynamics of this system, the r's represent trapping 
rates, e.g., r, on the rate at which electrons (n) are trapped by empty 
traps (t ) per electron per trap, and the g's represent release rates from 
the traps per trap. As we shall see, the key to determining the noise 
correctly within this model is to introduce source terms for each of 
these eight processes. 
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Fig. 2 — The eight trapping processes contributing to recombination. 

Even before we consider this problem quantitatively, we realize that 
correlations among fluctuations in the various particle densities will be 
complicated. For example, a positive fluctuation 8n in the number of 
mobile electrons will lead to positive fluctuations in the recombination 
directly via increased trapping in the t p traps and indirectly via in- 
creased trapping in the t„ traps, and the increased number of U traps. 
Using our microscopic approach we can deal with these correlations in a 
routine manner. 

The step-by-step application of our method to this problem is carried 
out in some detail in Appendix D. If we make use of the definitions 
Sn = 8?ii, 8p = 8n 2 , Bn tn = 8h 3 , 8n tp = 8ji 4 , Sn = (8>ii, 5?? 2 , 8?i 3 , 5n 4 ), 
then the cross-spectral density of 8n i} 8tij is given by 

sfrif) = £ tf*(«)fl*COJfS(«). (53) 

kl 

according to (124a). In (53), S d is the (4 X 4) cross-spectral density of 
the driving sources given by (124b) and M is the (4 X 4) matrix 
expressing the induced fluctuations 5n in terms of the spontaneous 
ones d : 

5n(a>) = M(o»)d(w), (54) 

according to (123). [Note from Appendix D that the vector d = (rfi, 
di, ds, dt) consists of four linear combinations of the eight fundamental 
noise driving terms that arise using this method. That these terms are 
mutually correlated can be seen from (124b). J Thus, from (53) we can 
calculate the spectral densities of interest. These include Su for the 
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spectral density of 5n, S 22 for that of 8p, and S i2 and £ 2 i for the cross- 
spectral densities between 8n and 8p, all of which contribute to the 
current fluctuations observed at the contacts of the device. Even more 
components of S 6a are needed to calculate fluctuations in recombination 
radiation, if such results are desired. 

We should stress that (53) contains all correlation effects (through 
S d ) and all relaxation effects (through M and M*) contained in this 
problem. If certain trapping or release rates are small relative to others, 
the contribution of such processes to S d can be neglected ; if M can be 
characterized by a single time constant t, MM* becomes proportional 
to (w 2 + 1/t 2 ) -1 , as is usually assumed. 8 The point to be made is that 
with no additional effort of a statistical nature we can solve device- 
noise problems which contain rather complicated statistical correla- 
tions. The key is to put the correlations in the dynamics of the problem 
of interest, thereby keeping the statistics simple. 

Finally, the reader is cautioned against making a normal mode 
analysis of the linearized noise equations and then introducing a 
noise-source term for each normal mode, since such terms will not be 
statistically independent in general. We note in concluding that, in 
obtaining (53), nowhere did we have to make use of normal-mode 
analysis. 

V. CONCLUSION 

In this paper we have discussed a straightforward, microscopic ap- 
proach which can be used to calculate the statistical fluctuations ac- 
companying any transfer process. Since nearly all charge-carrier ve- 
locity fluctuations can be characterized in terms of transfer processes, 
we, therefore, are able to calculate the statistics of the velocity fluctua- 
tions from those of the transfer fluctuations. Knowledge of the velocity 
fluctuations is often all the microscopic information that is needed to 
insert into the Langevin method or the impedance-field method to 
calculate device noise. (These methods explain in great detail how to 
convert velocity fluctuations into observed current and voltage fluctua- 
tions.) In so doing, we are able to insure that all important correlation 
effects and relaxation effects are included in the results. We further 
simplify the statistical portion of the calculation by separating the 
spontaneous from the induced fluctuations, and expressing the statistics 
of the latter in terms of the former. In this way, we have to calculate 
the statistics of only the uncorrected, spontaneous fluctuations from 
probability theory, while the more complicated statistics of the induced 
fluctuations can be obtained directly from those of the spontaneous 
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type. The correlation effects are thus shifted from the statistical portion 
of the problem to the dynamical portion. We have included (in Ap- 
pendix C) a derivation of an important result for the diffusion noise of 
"hot" charge carriers, that is, for charge carriers whose mobilities are 
field dependent. As the method developed here is readily generalized 
to nonstationary noise, our results will bo used extensively in treating 
noise in charge-transfer devices. 
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APPENDIX A 

In this appendix, we consider in detail a single elementary random 
process and calculate its statistical distribution. As we have shown in 
the text, many complicated random processes can be decomposed into 
simpler processes of the type considered here. Once the statistical 
distribution of these simpler processes is understood, one can determine 
the distribution of the more complicated process of interest. In the 
text, we have focused attention on autocorrelation functions, and on 
how to distinguish the autocorrelation function of a complicated 
process from those of the simpler processes of which it is composed. 
In Appendix B, we show how to obtain distribution functionals of 
complicated processes from those of simpler ones. 

The random variable whoso statistics we seek is 

d{t) = ZHt-U) - R"(t), (55) 

i 

where R"{t) (a specific function of time) is the mean rate of occurrence 
of the ti at time t. In other words, the U are independent random vari- 
ables, each specifying the time at which an independent random event 
occurs. To be a meaningful construct, it is necessary that R°(t) satisfy 

*°w»'^™' (56) 

that is, that the characteristic time in which R°(t) changes is much 
longer than the average time necessary for an event to occur. If (56) is 
satisfied, many events will almost surely occur in time intervals during 
which R"(t) changes only by a small amount. During these intervals, 
the statistics of the ti will be Poisson with mean rate R°(t). 

In probability theory, one often works with distribution functions. 
Thus, if a random variable x has a probability of p x (xi)dxi at x\, that 
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x is in dx h then the distribution function q p (k) of p x (x) is denned by 
q p (k) = / dxe ikx p x (x). (57) 

J— oo 

Since 

_^- ikx ^{x), (58) 

it is clear that q P {k) contains as much information about x as does 
Px(x). What is important is that it contains this information in a more 
convenient form, since, for example, all moments of x can be obtained 
at once from q P (k) by suitable differentiations with respect to A;. 

Here, we shall calculate the distribution functional 15 Q<i(k(t)) of 
d(t) defined by 

Q d [k(t)3 = /exp (i J'' dtk(t)d(t)\\ , (59) 

where the brackets denote averaging over the probability distribution 
of d(t). Rather than use the form of d(t) given in (55), let us use instead 
the form 

d(t) = E /(* - U) - R°(t), (60) 

where f(t) is a function of t satisfying 

/_" dtf{t) = L (61) 

Equation (59) then becomes 

Qd[k{t)-] = /exp U J' 2 dtk(t) (l)/(* - U) - R°(t)^]y (62) 

Let us now evaluate this average. 

To determine Qd[.k(t)^\ we note the following. 

(i) R°{t) is a specific function of time and, hence, 

exp (-i f dtk(t)R°(t)j 

can be factored out of the brackets. 
(ii) The U are independent and, hence, the average of the exponen- 
tial in J^,if(t — ti) can be factored into products of averages of 
an exponential in each f(t — U), each of the n averages being 
equal to all the others. 
(Hi) The probability that n events occur in the interval t\ < t < t% 
is Poisson. 
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From these three considerations, it follows that 
Qd[k(t)l = exp(-i J'* dtk(t)R°(t)\ 

X f Q P H [/exp (i j* dtk{t)f{t - *,)))]" , (63) 

where P„ is given by 

P„ = <r A A*/n! (64) 

and A is given by 

A - f' 2 R°(t)dt. (65) 

To evaluate the average in (63), we note that the probability that U 
occurs in dt' at I' is R°(t')dt'/\. Thus, 

/exp(t j h dt/c(t)f(t - ti)\\ 

= | f ' d«i exp (i j tl dtk(t)f(t - ti)\ R°(t>). (66) 

Inserting this into (63) and carrying out the sum on n, one obtains 
finally for Q d \Jc(t)2 

QiLk(t)l = exp I j\lt'R°{t') [exp (i J\ltk(t)f(t - *')) - l] J 

X exp (-i f h k(t)R°(t)dt\ (67) 

We note from (59) that if k(t) = 0, then Q d should equal 1. Since 
(67) satisfies this condition, we are assured that Q&[k{1)~] is properly 
normalized. As we shall see in Appendix B, from (67) we can calculate 
the distribution functional of nearly any function of d(t). 

It often happens that the particular process of interest involves 
numerous events, each of which contributes only a very small portion 
to the total fluctuation. In this case, we can expand the exponential in 
f(t — t') to quadratic order, obtaining for Qd[k(t)2 

<M>(0] - exp [- \ J' 2 dt'R°(t') (f*'dtf(t - «0*(0)*1 > (68) 

where we have assumed that f(t) is of sufficiently short duration that 

['* dt'R<>(t')f(t - t') - R°{t). (69) 
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Recalling finally that f(t - V) = 8(t - t') for the d{t) of (55), it follows 
that 

<M>(0] = exp (- \ J'* dtk(t)k(t)R°(t)\ (70) 

In the text, we made repeated use of the autocorrelation function of 
d, (d(h)d(t 2 )). This we may calculate from (67) with f(t) = 8(t), or 
from (70), the results being the same. To do this, we note that from (59) 



(d(h)d(t 2 )) = - 4- 



~ I QdtyiKt - h) + 728(1 - «,)]. (71) 

7i=0 072 | 7J =0 



Thus, from (70) we readily find that 

(d{h)d(t 2 )) = R°(h)b{h - U). (72) 

Although (67) is our most general expression for the statistics of d(t), 
(72) is usually all that is needed in noise calculations for devices. In 
simplest terms, it is the autocorrelation of the (shot) noise associated 
with the process of which d(t) represents the source term and R°(t) 
the average rate at time t. It is the building block from which most 
device noise can be constructed. For stationary processes, R°(ti) = R°, 
a constant depending on the noiseless solution. The spectral density 
of d under such conditions is a useful concept and is given by 



S d (f) = 2R°. (73) 



APPENDIX B 



In this appendix, we show how to obtain the distribution functional 
of a statistical process composed of a number of simple independent 
processes of the type discussed in Appendix A. Suppose that the process 
of interest, Bn(t) = n{t) — n°(t), is a linear functional of a number of 
d(t), \di(t), 1 ^ i ^ m\. Then 

8n(t) = JF n [di(0, • • •, d m (t), Q. (74) 

Since 8n is linear in the d„ we can write $„ in the form 

ff B = E Zm[.di(t), Q. (75) 

t-i 

The distribution functional Q,,[_k(t)2 of bn(t) is then given by 
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(76) 
(77) 

(78) 



(79) 



Q.[*(0] = (exp (i j\{t)Z n dt\\ 

= ft (exp (i J' 2 k(t)5nldi(t), Qdtj\ 

m 

where (77) follows from the independence of the di. 
A typical functional dependence of 8n(l) on d,-(Z) is 

8n(t) = £ a, [' dt'expf- J' dt"/ Ti (t")\ d&'). 

It follows from (78) that Q»<[&(0] is given by 
&.<[*«)] 
= /exp [» T' dflfc(*) /' dt'aMt') exp (- f <tt"/ *•*(*"))]) (80) 

= /exp Ti I' 1 dtaMt) f t2 dt'k(t') exp (- f d*"/T,-(o)"|\ (81) 

= Qn [a, j^ dt'k{t') exp (- J'' d*"/rf(*")]Y (82) 

If we use (70) for Q d „ then we find that 

Q ni [k{t)-] = exp j~ - \ J'* dt J* dtfk(t)k(t')Fi(t, *')] , (83) 

where 

JtfJ, *) = a 2 <exp [- f* m "dt"/Ti(t")] 

/Train I /• T m |n I 

dr«?(T) exp -2 J dt"/ Ti (t") (84) 

and where r max = max (2, t'), r min = min(<, t'). Inserting (83) into (78) 
determines Q n as desired. 

From Q n we can determine the autocorrelation function of n(t). 
First, we note from (74) and (76) that 



(8n(h)8n(t 2 )) 



_a_ 



d 

7,-0 ^72 



72=0 



QLyi8(t-h) + y£(t - fa)]. (85) 
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Using (78) and (83) this becomes 

(Sn(ti)8n(t2)) = Z Fitth fa), (86) 

» 

where F t is given in (84). 

If r, and R"i are independent of t, an interesting result is obtained. 

By performing the integrations in (84), it follows at once that 

Fi(h, l 2 )/a 2 i = (RlTi/2) exp (- \t 2 - h\/n), (87) 

where the correlation of the fluctuations contributing to 8n is readily 
apparent. Expression (87) occurs often in the theory of stationary 
processes. [Note in (87) that in the limit n — » 0, F"{t\, ^)/a 2 i approach 
Rtfiti - h) as expected from (72).] 

APPENDIX C 
Fluctuation-Dissipation Theorem 

Relations between different physical phenomena have always at- 
tracted considerable interest, and justifiably so. No exception is the 
Einstein relation between the mobility n and the diffusion coefficient D : 

„-g- (88) 

In this expression, m relates the carrier's velocity v to the electric field 
E acting on the carrier according to 

v = nE, (89) 

and D relates the mean-square distance which the carrier diffuses in 
equilibrium to the time t in which it has been diffusing according to 

<(*, - x )*) = 2Dt. (90) 

Although usually not emphasized, all derivations of (88) are based on 
small fluctuations from equilibrium (in the absence of driving forces). 
In the presence of driving forces, one must be careful to use the ap- 
propriate small-signal quantities when relating diffusion to transport. 
One of the purposes of this appendix is to explain how this can be done 
in some cases. 

It is very important at the outset to understand the physical origin 
of the Einstein relation, as well as that of the more general fluctuation- 
dissipation theorem 16-18 (FD theorem). This is especially important 
for noise theory since it often enables one to express fluctuation proper- 
ties in terms of transport properties. This is a valuable aid since trans- 
port properties have already been carefully studied in attempts to 
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understand the noiseless behavior of the device. In addition, an under- 
standing of how the microscopic noiseless motion of a single carrier 
contributes to the noiseless device current can be carried over at once 
to calculating how the microscopic fluctuations contribute to the device 
noise. Thus, we ask the question, Why are fluctuations and dissipation 
so closely related? 

What is the physical origin of the fluctuations undergone by charge 
carriers ? The answer clearly lies in the scattering of all sorts that such 
carriers experience within the material. What is the physical origin of 
the damping force experienced by charge carriers? The answer clearly 
is the same. This connection may be phrased in several ways. A fluctua- 
tion corresponds to the response of the carrier to a random force. The 
velocity-field dependence is a similar relation of response to applied 
force. Alternatively, under steady-state conditions, the average gain 
in the energy of a carrier from the material due to fluctuations must be 
dynamically balanced by the loss of energy due to damping. The gain 
and loss are, therefore, closely linked. Still another way of seeing the 
connection is to note that a fluctuation is a departure from equilibrium, 
such as when a small, disturbing probe force is applied to the system. 
However it is viewed, a close connection between fluctuation and 
dissipation must clearly exist, which we shall derive below. 

In the impedance-field 4 method, we arrive at the following equation 
for the voltage spectral density of the device in terms of the elementary 
thermal- velocity fluctuations of the carriers : 

S(BV N) f) = £ f ' dx\V 1 Zwz\*4#n(x)D t (to it f), (91) 

1=1 J 

where q is the elementary carrier charge, n(x) is the carrier density, 
Znx is the impedance field between the contact at N and the field point 
at x, and 

S.(x, /) = 4D x (5y„ /) = 4Re /"" e iut (8Vi(t)Bvi(0))dt (92) 

Jo 

is the velocity spectral density in the ith. ( = x, y, z) direction. [See 
eqs. (54), (56), (27), (28), and' (35) of Ref. 4.] The quantity D x (8 Vi , f) 
is also referred to as the diffusion of 5v, at frequency / in the region of 
x for a single carrier. In discussing the impedance-field method for the 
case of v = y.E, where y. is a constant, it is found that 

D x (h>*,f) = (y)" < 93 > 
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for / much less than the scattering frequency. In what follows, we use 
the FD theorem to prove that if n = n{E) = v(E)/E and if an effective 
temperature T can be defined, then 11 

[Intuitively, (94) follows from (93) if we recall that fluctuations in 
velocity result from spontaneous fluctuations in the electric field, and if 
v = v{E), then 8v = [dv(E)/dE']-8E = (d[fi(E)E3/dE)'8E.l In the 
Langevin' 1 method, diffusion noise enters through a spectral density of 
the form 

S h (x, x', /) = 4qhi(x)D n 8(x - x'), (95) 

where again D„ is the diffusion constant, as defined in (42) for fre- 
quencies much less than the scattering frequency. Noting how (42) 
is related to (8v(ti)8v(t 2 )), we note that, once we have shown that (94) 
is true, it will follow at once that 

S h (x,x', /) = 4 q hi(x) y ^ 8(x - x'). (96) 

The first quantum-mechanical derivation of the FD theorem is 
usually attributed to Callen and Welton. 10 This subject has subse- 
quently been treated in greater detail in Refs. 19 and 20. For the sake 
of completeness, we shall redcrive the FD theorem here primarily for 
the purpose of calling attention to its application to cases of steady- 
state but nonequilibrium conditions. We conclude this appendix with 
the derivation of more general relations for fluctuations and dissipation 
valid for any energy distribution (especially nonthermal) of the states 
of the system. These relations show how closely the two are related 
even under nonthermal conditions. 

Let us suppose that we have some system which can be described 
by a Hamiltonian H*, which includes applied electric and magnetic 
fields giving rise to currents, etc. Let f t be an operator whose expecta- 
tion value we seek as a function of time in response to a unit impulse 
in a probe force /,-, which enters the total Hamiltonian // according to 

H = H°- fjfj, (97) 

where r, is another operator. Then the FD theorem states that 

B.pr.MJ-lmg^a ^ + gjlg , (98a) 
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where /3 = 1/kT, T is the effective temperature of the system, if such 
can be denned, 

KM = I* dtkiiiW, (98b) 

Jo 

MO = *<C><(0, ^(0)]+) - KM04CO) + *(0)MO>, (98c) 

Xij(w) = [ X dtxa(t)e iat , (98d) 

and 3^(0 is the impulse response of r, to the unit impulse /,■ in (97). 
From (98b, c) it is clear that 4 Re [Ka( a )l = 5 ,->(/), the cross-spectral 
density between r,- and r>, whereas from (98d) it follows that 
Im £Xij(u)2 is the absorption through r,- of the energy put into the 
system via the coupling (interaction energy) — /,-ry. 

Returning for the moment to (92), we note that if we let r, = fj 
= 8Vi = 8x (for i = x), in the classical limit (h — »0), we obtain from 
(98a) 

S v (x, f) = 41m [X«(»)]fc!T/ft>. (99) 

After we derive (98a), we shall show that for w much less than the 
scattering frequency 

Xii (w) = iunac/q, (100) 

where /Xoc = d\_n{E)E~]/dE. The important result (94) then follows at 
once. 

We shall derive (98a) by evaluating A",/(a>) and Xij{u>) and com- 
paring the results. By definition 

MO = \ (ri(t)n(0) - r>(0)r,(0), (101) 

where we assume the impulse occurs at t = 0. Clearly Xa(t) = 0, t < 0. 
Thus, 

Xa(fi») = \ r die- 1 £ (l|Mm)<m | >/ 1 1> (r**« - «"**■) 

x e-i/*{**-J!i)«/£ c -/»*i, (102) 

which follows if we assume that the system described by #' can be 
characterized by a temperature T. This insures that the probability 
that the system is in eigenstate \f/„ with eigenenergy E n is 
oxp { — (3E„)/N, where N is the normalization factor used in (102). 
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Performing the integration on t, we obtain 

lm 

X (-l)O - (E m - E t ) + Kj-ytf, (103) 
and, finally, taking the imaginary part we find that 
lm [?«(«)] 
= 7T Z (**■• - «"*■-) <Z 1 r,- 1 wi><m | fj \ l)6[faa - (E m - E,) ~\/N. (104) 

lm 

If now we compare the definition of fcy(<) (98c) with (101), it follows 
at once from making the appropriate changes in (103) that 

KiM = is (e-^ £ ' + *-«-)<! |Mm)(m|MD 

X (-1)[> - (*« - #i) + uQr l /N- (105) 
And upon taking the real part of (105), 

Re [*<,(•)] = * \ 2 (e-^' + *-**-) 

X <i|/ f |m)<m|/y|Z)8C^ - (£ m - #,)]/# (106) 

is obtained. In (104) and (106), we may replace E m in the exponent by 
(Ei + ^w) owing to the presence of the delta function. This then 
permits the factor [1 — exp (— /3/zaj)] to be pulled out of the sum in 
(104) and [1 + exp (-jStfw)] to be pulled out in (106). The FD 
theorem (98a) follows at once. 

To derive (100), we must look more closely at the exact response 
x\j of fi to a general time dependent /,- in H = H s — fjfj. If H is the 
Hamiltonian, then 

xh(t) = (fih = Tr [MOU/Tr [p(01 (107) 

where p(t) is the density matrix of the system at time t. The density 
matrix p(t) at t can be obtained from that at t\ according to 

p{t) = P&tMtJPt&h), (108) 

where the propagator P{t, t\) is defined by 

P(t,h) = exp(-!yVff r ), (109) 

and where we are using the Feynman 21 ordered-operator notation. [li 
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h is prior to the turning on of / and the system can be described by an 
effective temperature T, then p(h) = exp (— {$H").'] 

The specific problem in which we are most interested is the case 
where f i} fj represent velocity fluctuations from the expectation value 
of the velocity v of a carrier under the influence of an electric field of 
arbitrary strength contained in H". Let us place ourselves in a reference 
frame drifting with velocity v . In such a frame, we may calculate 
either the velocity fluctuations fc»(t) from the noiseless motion or the 
response Xu(t) to a probe force /,. If we were applying a small, ac 
electric field e(t) to the carrier, then H = H" — qe(t)£, where A is the 
position operator corresponding to the x coordinate (in the drifting 
frame). Also, we know that the velocity response to such a field is 
given by 

8v(t) = ldv(E)/dE~\e{t) (110) 

or, taking Fourier transforms, 

*»(») = \dbi(E)E]/dE}e(u>), (111) 

since v(E) = n(E)E, n being the mobility. However, X x i(w) arises 
from the velocity response to H = H" — f(t)x. We can relate f(t) to 
e(t) if we note that from (109), f(t) enters in the form 

exp (- \ f dT/(r)i) = exp (- \ f' *r(-)/(r)s) , (112) 

where we can take fit) = for simplicity. Thus, qe(t) = —fit), or 
qe(co) = ioif(co). Also, using (111), we obtain 

Y \ = gp (") = *"(") ^ - d\ji(E)E] ia) 
***{<*) - y (w) ; = e(w) q - dE q 

= iufiac/q, (113) 

which proves (86). Our principal result (94) is therefore demonstrated. 

We promised that we would conclude with more general relations for 
fluctuations and dissipation valid for any probability distribution of 
the eigenstates of a system according to eigenenergies. We proceed as 
follows. 

Let the (normalized) probability that the system be in a state of 
eigenenergy E be f(E). Choose the scale of energy such that the lowest 
eigenenergy is zero (0). Take the Laplace transform of f(E) to obtain 

F(s) = r dEf(E)e- E (114a) 

Jo 
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and, of course, 

1 /■o+t'oo 

f(E) = —-. J dse' E F(s). (114b) 

[For a thermal distribution, f(E) = exp (—0E)/N and F(s) = 
(s + /3)-yiV.] It follows from (104) and (106) that 

| Re [£«(«)] = 2^- / _ . <«?(«, «) (1 + «•*") (H5a) 

and 

1 ra+ito 

Im [Z (o))] = ~ / dsG(s, «)(1 - 6'""), (115b) 



where 
G(s, u) = wF(s) E e«*'<J|*|m)<w|^|J)a[>> - (#m - Bi)J (115c) 

The difference between (115a) and (115b) consists of only one sign: 
plus in (115a), minus in (115b). For the special case of a thermal 
distribution, this sign difference yields the ratio factor that appears in 
(98a). In general, while such a simple relation between Ky and Xa is 
no longer valid, it is clear from (115) that fluctuations and dissipation 
have a common origin. These are the more general relations we 
promised. 

APPENDIX D 

The purpose of this appendix is merely to carry out the routine 
mathematical steps necessary to arrive at the spectral densities of 
recombination-generation noise. These results are referred to in Section 
IV. 

We proceed as follows. We can abbreviate our exposition since we 
have been through the necessary steps several times in the previous 
sections. The steps are as follows: 

(i) Write the microscopic kinematic equations : 

dn 

-r t ■ = - E «P - t ten (m)2 - £ ap - t tpn (m)2 

at tn m 

+ E S[t - tntn(m)2 + E 8ft - tnt (m)2, (116a) 

m m 

% = - E ap - t t0P {m)-\ - e i[i - *!**(»»)] 

+ £ 8p - Wro)] + I«[f- «p to (m)l (116b) 

?n m 
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^r = + E i[« - * to »(m)] - E«p - *n«»(w)] 

- E «[< - <£»p(w)] + E «B - ^ to (m)l (H6c) 

Hi rn m 

+ L5[/- t lop (m)2 - E «D - W™)1 (116d) 



dn« _ rf/iin dn t p 
dt dt dt 

(ii) Rewrite in terms of sources and responses : 



(116e) 



dn 

= —r, on nt n — r lpn n tp n + g 

- d t0 n — d tp „ + d nln + d„, , (117a) 



-jj = —r, on ?i to n — r lpn n lp n + g n tnn tn + g nto ?ito 



dp 

jT — TtopfltoP TtnpTlinP "T Qplpf^lji T" Q P lo1llo 

— diop — dmp + dptp + d pt0 , (117b) 
-g- =r t0 nn to n — g n i n n tn — r t n V n tn p + g P ion lo 



-j~ = — r tp „n tp n + g„ lo n to + r lop n lo p — g ptp n tp 



drito _ _ <in <n _ drii p 
dt dt dt 



+ dton — dntn ~ d tn p + d pto , (117c) 

pUtoP — gpipUtp 

— dtpn + d nt o + dio P — dptp, (117d) 

(117e) 



where a typical d function is defined as 

d l0 n = E *(' — t'on) — r,o n n,„n (118a) 

m 

and is equivalent for small fluctuations to 

d lon - E *(< - C) - O&n* (118b) 

(m) Write each variable as the sum of a noiseless contribution and a 
noise contribution, and linearize (117) to obtain the following nonlinear 
equations for the noiseless solution : 

n° = —rton>Cn° — r lpu n ipn + gnt n n°, n + g n ton? , (119a) 

p° = -r,o P nt°op - r„ ip n° m p + g plp nt p + g pto n° , (119b) 
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n°tn = r tannin - g n tnn tn - r tnp n%p -f g p i n%, (119c) 

n% = -r tpn n lp n + gnt o n° t0 + rt op n lo p° — 9vt P n° P) (119d) 

nl = -n?„ - n° lp , (119e) 

and to obtain the following linear equations for the noise solution : 

Sn = —r t on(n?o&n + 8n to n°) — r tP n{n tp bn + 8n lp n°) 

+ g n tn8n t n + gmo8n, — cZ, „ — di pn + d„,„ + rf n /o, (120a) 

5p = -r top (n%8p + 5n, p°) - r, np (n?„5p + 5n (n p°) 

+ g P t P 8n lp + g p ioSn to — rf( p — d<„ P + dpfp + d pto , (120b) 

S?i<„ = r, on (n? 5n + 5n< n°) — g nt n8n t n — r tnp (n? n 8p + 8n tn p°) 

+ g P to8?llo + dfon — dntn — d*np + d p «,, (120c) 

5ri ip = —r tpn (n1 P 8n + 6n, p n ) + g n to8n to + r, op (n? 8p + 5n to p°) 

— g p t P 8n lp — d tp n + d„ to + diop — d P < P , (120d) 

fin< = — 5rij„ — 5n ep . (120e) 

(w) Finally, solve (120) for 8n, 8p, 8n tn , 8n tpi and 8n l0 in terms of 
dton, • • • . This latter step involves solving five equations for five un- 
knowns. If we are interested in the stationary solution to (119) 
(n° = p° = n° ln = n° p = 0), then the coefficients of the 8n, 8p, 
(120) are independent of time and a Fourier analysis of (120) is most 
expedient. Using (120e) to eliminate 8n l0 , one obtains the following set 
of equations : 

iw8n = —R„8n + R n tn8n tn — Rnt P 8n tp + d\, (121a) 

iu8p = —R p 8p — R ptn 8n t n + R ptp 8n tp + d 2 , (121b) 

iw8nt„ = —Rtn8n ln + Rm n 8n — R tnp 8p — R ln t p 8n tp -f- d3, (121c) 

io>8n tp = —Ri p 8n tp — R tpn 8n + R tpp 8p — Rt pt n8n tn + c^, (121d) 

where 

dl = ~dton — dtpn + dnm + dnio, (122a) 

d 2 = — d (op — d, np + d p t p + dp< , (122b) 

^3 — dton — dnm — d lnp + d pt0 , (122c) 

and 

di = —d tnp + dn«o + d<op — dp, p . (122d) 

By combining the coefficients of the 8n, 8p, • • • in (120), the reader can 
determine the R n , R P , • • • rates in (121). For example, R n = r l0 nn°o 
+ rt pn n° t p, etc. The important point to be made here is that the d\, dz, 
d 3 , di are correlated, that is, they are not mutually independent. Thus, 
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one cannot a priori introduce independent d h d 2 , d 3 , d. t into equations 
of the form of (121) at this stage of the calculation. Using our micro- 
scopic approach, however, these source terms and their correlations 
arise in a well-defined manner. 

If we let 8n = 8n h 8p = 8n 2 , 8n, n = 8n 3 , 5n tp = 5w 4 , Sn = (8n h 8n 2 , 
8n 3 , 871a), d = (di, d 2 , d 3 , dj, then (121) can be written in the form 

iultak = R5n + d (123a) 

and solved for 5n at once : 

8n(w) = M(w)d(«), (123b) 

where 

M(w) = (iwl - R)" 1 . (123c) 

The cross-spectral density of 5n„ 6nj is then given by (« = 2ir/) : 

89(f) = L MMSUflMfa), (124a) 

kl 

where £&(/) is the cross-spectral density of d k , dr. 



S" = 



•J|oiiTU|piiTWw»T'Jni 



Sd Qd 
ton &ntn 

l_ *Spn I &nto 

Sd Qd 

ton &nln 

OmprtJ Diu 

Oton ~T&iUn T ®tnv i&pto 







&lnp TOinpT Optp ~T »J pin 

Olnp\&pto 

Sd Qd 

lop &plp 

Sd I Qd 
lpn\&nto 

Sd Qd 

top &ptp 



Stp r -\-£>nto rOfop + Of 



(124b) 



and 



&ton — <uf to tin", &top — £T topP , &tpn — *Ttpn"i j 

Sfnp = 2rtn P p , Sim = 2g n in, S% p = 2g P t P , S nl0 = 2g nt0 , 

&plo = ^Qpto 

from (73). Using (123b, c) and (124a, b), the cross-spectral density of 
8ni, 8iij including all correlation effects can be readily calculated. Of 
greatest interest are the spectral densities of 8n and 8p given by Silif) 
S'£(f), and of 8?i, 8p given by Si§(/) = Sfi(/). These are the terms in- 
volving mobile carriers, which contribute to the current fluctuations 
seen at the contacts of the device. It should be clear that trapping can 
be treated by using an approach similar to the above. 

There is little point in carrying out the algebra implicit in (124a). 
Evaluating the inverse of the 4 X 4 matrix in (123c) and carrying out 
the sum in (124a) can be done readily, although it is somewhat tedious. 
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What we have accomplished is to express the cross-spectral densities of 
the electrons, holes, and various trapping states in terms of known, 
easily calculated, cross-spectral densities of the source-driving terms. 
And in so doing we have included all important correlation effects in a 
very natural way. 
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